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Abstract: 

I calculate the form factors describing semileptonic and penguin-induced decays of B mesons 
into light pseudoscalar mesons. The form factors are calculated from QCD sum rules on 
the light-cone including contributions up to twist 4, radiative corrections to the leading 
twist contribution and SU(3)-breaking effects. The theoretical uncertainty is estimated to 
be ~ 15%. The heavy-quark-limit relations between semileptonic and penguin form factors 
are found to be valid in the full accessible range of momentum transfer. 



^Heisenberg Fellow 



1. Decays of B mesons into light mesons offer the possibihty to access the less well known 
entries in the CKM quark mixing matrix such as Vub and Vts- The measurement of rare 
penguin-induced B decays may also give hints at new physics in the form of loop-induced 
effects. With new data of hitherto unknown precision from the new experimental facilities 
BaBar at SLAC and Belle at KEK expected to be available in the near future, the demands 
for accuracy of theoretical predictions are ever increasing. The central problem of all such 
predictions, our failure to solve non-perturbative QCD, is well known and so far prevents 
a rigorous calculation of form factors from first principles. Theorists thus concentrate on 
providing various approximations. The maybe most prominent of these, simulations of QCD 
on the lattice, have experienced considerable progress over recent years; the current status 
for B decays is summarized in [Q]. It seems, however, unlikely that lattice calculations will 
soon overcome their main restriction in describing b —>■ u and b s transitions, namely 
the effective upper cut-off that the finite lattice spacing imposes on the momentum of the 
final-state meson. The cut-off restricts lattice predictions of B decay form factors to rather 
large momentum transfer of about 15 GeV^ or larger. The physical range in B decays, 
however, extends from to about 20 GeV^, depending on the process; for radiative decays 
like B K*'-f it is exactly GeV^. Still, one may hope to extract from the lattice data some 
information on form factors in the full physical range, as their behaviour at large restricts 
the shape at small via the analytical properties of a properly chosen vacuum correlation 
function. The latter function, however, also contains poles and multi-particle cuts whose 
exact behaviour is not known, which limits the accuracy of bounds obtained from such 
unitarity constraints and until now has restricted their application to B ^ n transitions 
0, ^ . The most optimistic overall theoretical uncertainty one may hope to obtain from this 
method is the one induced by the input lattice results at large g^, which to date is around 
(15-20)% 0, 0. A more model-dependent extension of the lattice form factors into the low 
region is discussed in 0. 

An alternative approach to heavy-to-light transitions is offered by QCD sum rules on 
the light-cone. In contrast and complementary to lattice simulations, it is just the fact that 
the final-state meson does have large energy and momentum of order ~ mBf^ in a large 
portion of phase-space that is used as starting point (which restricts the method to not 
too large momentum transfer, to be quantified below). The key-idea is to consider b ^ u 
and b s transitions as hard exclusive QCD processes and to combine the well-developed 
description of such processes in terms of perturbative amplitudes and non-perturbative 
hadronic distribution amplitudes (see also [0 for a nice introduction) with the method 
of QCD sum rules [§] to describe the decaying hadron. The idea of such "light-cone sum 
rules" was first formulated and carried out in M in a different context for the process 



S p7, and its first application to B decays was given in [T^. Subsequently, light-cone 
sum rules were considered for many B decay processes, see [11, |l^ for reviews. As light- 
cone sum rules are based on the light-cone expansion of a correlation function, they can be 
systematically improved by including higher twist contributions and radiative corrections to 



perturbative amplitudes. The first calculation in |T0[ was done at tree-level and to leading 
twist 2 accuracy. In |]13|, twist 3 and 4 contributions were included, and in [|^, one- 
loop radiative corrections to the twist 2 contribution to the form factor /I were calculated. 
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In [TO], the corresponding radiative corrections to the decays of B mesons into the vector 



mesons p,K*,(j) were calculated. In [0], the scalar form factor was calculated at tree- 
level. In the present letter, I calculate the remaining radiative corrections to all semileptonic 
and penguin B —>■ ir and B K transitions and present new and more accurate results for 
the corresponding form factors. 

2. Let me begin by defining the relevant quantities. The semileptonic form factors are 
defined as {q = Pb — p) 

{P{p)\qi,b\B{ps)) = f^iq') i^{pB+p), -^^^^^q,] + ^^^^^ fo^i^') q„ (1) 

where P stands for the pseudoscalar meson ir or K and q = u for the tt and q = s for the 
K. The penguin form factor is defined as 

{P{p)\qa,,q''{l + ^,)h\B{pB)) ^ {P{p)\sa,,q%\B{pB)) 



(2) 



' itlb + rriK 

The physical range in is < < (m^ — mpY' . Although there are of course no semilep- 
tonic decays B — »• Kev, the above form factors contribute to, say, B —>■ Kii. Recalling 
the results of perturbative QCD for the vr electromagnetic form factor as summarized in 
0, one may suppose that the dominant contribution to the above form factors be the ex- 
change of a hard perturbative gluon between, for instance, the u quark and the antiquark; 



this possibility was advocated in [|1^. This is, however, not the case, and it was pointed 
out already in Ref. |TD[ that the dominant contribution comes from the so-called Feynman 
mechanism, where the quark created in the weak decay carries nearly all of the final-state 
meson's momentum, while all other quarks are soft, and which bears no perturbative sup- 
pression by factors a^/vr. In an expansion in the inverse b quark mass, the contribution 
from the Feynman mechanism is of the same order as the gluon-exchange contribution with 
momentum fraction of the quark of order 1 — AqcB/'^b, but it dies off in the strict limit 
— > oo due to Sudakov effects. This means that — unlike the case of the electromagnetic 
TT form factor — knowledge of the hadron distribution amplitudes 

~ / dk]_'^{u,k±), 
Jo 

where \Ef is the full Fock-state wave function of the B and vr(i^), respectively, u is the 
longitudinal momentum fraction carried by the {b or u{s)) quark, k± is the transverse quark 
momentum, is not sufficient to calculate the form factors in the form of overlap integrals 

F ~ ^ dudv4>l(^j^^{u)Ti,^rd{u,v;q^) 4>b{v) 

(with Thard oc ds)-^ Instead, in the method of hght-cone sum rules, only the light meson 
is described by distribution amplitudes. Logarithms in k± are taken into account by the 

^Note also that not much is known about 0b, whereas the analysis of light meson distribution amplitudes 
is facilitated by the fact that it can be organized in an expansion in conformal spin, much like the partial 
wave expansion of scattering amplitudes in quantum mechanics in rotational spin. 
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evolution of the distribution amplitudes under changes in scale, powers in k± are taken 
into account by higher twist distribution amplitudes. The B meson, on the other hand, 
is described, as in QCD sum rules, by the pseudoscalar current di'y^b in the unphysical 
region with virtuality — rnl ~ 0{mb), where it can be treated perturbatively. The real B 
meson, residing on the physical cut at = m%, is then traced by analytical continuation, 
supplemented by the standard QCD sum rule tools to enhance its contribution with respect 
to that of higher single- or multi-particle states coupling to the same current. 

The starting point for the calculation of the form factors in (P and (0) is thus the 
correlation functions {Jb = di'y^b): 



CFv 



(3) 



CFn 



(4) 



which are calculated in an expansion around the light-cone = 0. The expansion goes in 
inverse powers of the b quark virtuality, which, in order for the light-cone expansion to be 
applicable, must be of order mj,. This restricts the accessible range in to — ~ 0{mi,) 
parametrically. For physical B mesons, I choose ml — q^< 17GeV^. Note also that for very 
large the influence of the next nearby pole [B* for /^) becomes more prominent. 

It proves convenient to perform the calculation for an arbitrary weak vertex F = {7^, 
ct^uQ'^}, which, neglecting for the moment radiative corrections, yields: 



CFr = — f^du 
4 JO 



iu)TriTS,iQ)p) 




6 OQa 



dvg2{v, 
J dQaoQ 



-Tr{TS,{Q)p)-g2{u)^Tr{TS,{Qho 



+ 4 f'dv 



4 Jo Jo 



4/3. 



v{pq)(t)^M)Tr{Vp) + (20^ (a) - <P\\{a))Tr{V{q + mb)p) 



+ 2v [(P\\{a)Tr{Vpq) - 2(pg)0x(a)Tr(F)} + {20^ (a) - <P\\{a)]Tr{V{q + mi,)p) 



+ ^iv(t)^{a)Tr{Va^p)qyP 



Explicit expressions for n± and Ft were already obtained in 



nit 



mt 



■up%- 



uq"^, Va = daida2da36{l — ai — 02 — cts) and s 



(5) 

Here Q = q + up, 
ml-{q+{ai + 



va^)pY] Sb{Q) = {Q + mb)/{—s) is the b quark propagator. In the above expression, 0^ /^ 
is the leading twist 2 distribution amplitude, (pp and (pa are the two-particle distribution 
amplitudes of twist 3, gi and g2 those of twist 4, all of which are defined in [|19|. The twist 3 



3 



and 4 two-particle distribution amplitudes are determined completely in terms of the twist 3 
and 4 three-particle distribution amplitudes 03^, and 0|| j_ [|T9]. Note that in the above 
expression corrections in the light meson mass are neglected (m^/(m„ -|- m^), however, is 
expressed in terms of the quark condensate and taken into account). Their inclusion, of 
potential relevance in B ^ K transitions, is not straightforward and requires an extension 
of the method developed in Ref. [|l^ to include meson- and quark-mass corrections in the 
twist 4 distribution amplitudes. According to [^, the numerical impact on the form factors 
is small, around 5%, and most pronounced at large g^. 

3. It is convenient to calculate also the radiative corrections for arbitrary weak vertex. To 
twist 2 accuracy, the light quarks are massless and carry only longitudinal momentum. The 
one-loop calculation does not occasion any particular technical complications, but results 
in bulky expressions which I refrain from quoting here. The general structure is, as to 
be expected, similar to that for the form factor obtained in |T^. The separation of 
perturbative and non-perturbative contributions introduces an arbitrary logarithmic (infra- 
red) factorization scale. The condition that the correlation function be independent of 
that scale leads to an evolution equation for the distribution amplitude, which was first 
derived and solved in to leading logarithmic accuracy. In the present context, with full 
0{as) corrections to the perturbative part included, one has to use the next-to-leading order 



evolution of the distribution amplitude, which was derived in closed form in ||21]|. A natural 
choice for the factorization scale is the virtuality of the b quark, /if-^ ~ urrib. For technical 
reasons it is, however, more convenient to choose a fixed scale like fif^ = m'^^—ml, which is of 
the same order. The numerical impact of changing the scale is minimal.^ The penguin form 
factor depends also on an ultra-violet scale, the renormalization-scale of the local operator 
QO'fj.uQ'^b appearing in the effective weak Lagrangian. A natural choice for this ultra-violet 
scale is /luv = mb- 

As for the radiative corrections, it turns out that they are dominated by the correction 
to the pseudoscalar B vertex, which, as discussed below, yields large cancellations against 
the corresponding corrections to the leptonic B decay constant /b- 

4. Let me now derive the light-cone sum rules. The correlation functions CFp, calculated 
for unphysical p%, can also be written as dispersion relations over the physical cut. Singling 
out the contribution of the B meson, one has for instance for 11+: 

CFn, = — f+i'f ) — o T + higher poles and cuts, (6) 

where fs is the leptonic decay constant of the B meson, fB^% = ''^b(-B|j^|0). In order 
to enhance the ground-state B contribution to the right-hand side, I perform a Borel- 
transformation: 



^This is in contrast to the tt electromagnetic form factor, which is rather sensitive to the shape of the 
distribution ampUtude near the end-points. 
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with the Borel parameter M^. The next step is to invoke quark-hadron duahty to approxi- 
mate the contributions of hadrons other than the ground-state B meson, so that finally 



B CFn, = ^ ^ U{q') e"™-/*^^ + - Hds ImCFn,(^) expi^s/M'). (8) 

This equation is the light-cone sum rule for /+; those for /o and /t look similar. Here Sq 
is the so-called continuum threshold, which separates the ground-state from the continuum 
contributions; Sq and are in principle free parameters of the light-cone sum rules, but 
they can be fixed by requiring stability of the sum rule under their change. In the present 
context, one can decrease their infiuence considerably by also writing fs as a QCD sum rule, 
depending on the same parameters Sq and M^. From the analysis of the latter sum rule, one 
finds So ^ 34GeV2 and GeV^. The resulting value for fe is (150-200) MeV, in 

perfect agreement with the results from lattice simulations. This procedure makes the form 
factors largely independent of nib, so and M^; the remaining dependence will be included 
in the error estimate. Note also that subtraction of the continuum contribution from both 
sides of (ID introduces a lower limit of integration u > {ml — 5^)/(so — 9^) in (^, which 
behaves as 1 — Aqcd / for large nib and thus corresponds to the dynamical configuration 
of the Feynman mechanism. 

Let me now specify the non-perturbative input. For the b quark I use the one-loop 
pole mass = (4.8 ± 0.1) GeV, which is consistent with a recent determination from the 
T mesons p2|]. For the light mesons, the distribution amplitudes need to be specified. 



Fortunately, conformal symmetry of massless QCD combined with the nonlocal string op- 



erators technique developed in provides a very powerful tool to describe higher twist 
distribution amplitudes in a mutually consistent and most economic way (see for a 
detailed discussion). The determination of the relevant non-perturbative parameters from 



QCD sum rules was pioneered in [^. In |jT9l, the twist 3 and 4 vr distribution amplitudes 



were obtained including contributions up to conformal spin 11/2 in terms of 6 indepen- 
dent non-perturbative parameters whose values were determined from QCD sum rules. The 
leading twist 2 distribution amplitude, on the other hand, can be expanded in Gegenbauer 
polynomials Q : 



= 6m(1 -u)il + J2 a^ifi)Cf\2u - 1) . (9) 



■y3/2 

\ i=l 

The Gegenbauer moments Oj renormalize multiplicatively. For vr, all odd moments vanish 
because of the tt's definite G-parity. In practice, one truncates the expansion after the first 
few terms, 

07r,ii-(n) = 6u(l -u) + ai{ij)Cf'^{2u - 1) j , 

with n = 4 (for vr) or n = 2 (for K) . I will discuss the impact of this truncation on the form 
factors later. For now, I use the vr distribution amplitude as obtained in p6| (see also p7[|), 

alii GeV) = 0.44, a^(l GeV) = 0.25. (10) 
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For the K, on the other hand, the non-zero value of the strange quark mass induces non- 
vanishing values of the odd moments. I use 



af (IGeV) = 0.17, (1 GeV) = 0.2, (11) 



where the first value was obtained in ||2^ and the second one comes from an analysis of the 



sum rule for the vr in due account being taken of SU(3)-breaking effects. 



The results are displayed in Fig. 1. The form factor coincides with the one obtained in 
T5[| . I plot each form factor using the twist 2 distribution amplitudes as specified above and 
with and without 0{as) corrections, and also using the asymptotic distribution amplitudes 
(pTT,{o) and (pK,{i) to illustrate the impact of non-asymptotic contributions. The plotted curves 
were obtained with nib = 4.8 GeV, Sq = 33.5 GeV^ and = 6 GeV^. The distribution 
amplitudes are evaluated at the scale //^ = — ^1 —'■ lA- Apparently, the net effect of 
radiative corrections on the form factors is rather small. This is due to an effect already 
observed in |T^: the radiative corrections to the QCD sum rule for are rather large, 



which is due to the large vertex corrections to the pseudoscalar B vertex. In the radiative 
corrections to the light-cone sum rules, the same vertex appears with corrections of similar 
size, so that they cancel between left- and right-hand side of (||), leaving a net effect of 
around 10%. 

For all form factors, the effect of three-particle twist 3 and 4 quark-gluon contributions 
(and their induced effects in the two-particle distribution amplitudes) are small (~ 5%), so 
that the considerable theoretical uncertainty of these terms does not play. This also shows 
that the expansion in contributions of increasing twist is under good control. The remaining 
twist 3 contributions are proportional to the quark condensate, which, as already noted in 



T5[| , introduces only a small uncertainty in the final results. 

As is expected from the definition of /o, which refers to a scalar current, it increases less 
sharply in than the other form factors. A good parametrization for the dependence 
can be given in terms of three parameters as 

F{(f) = ^ . (12) 



\ — ap {q^ /m\) + bp {q^ /m 



B) 



The parameters are given in Table 1 for central values of the input parameters. For com- 
parison, I also give the results for quoted in [jl2l and obtained in JTsf , the latter being 



obtained in leading-logarithmic accuracy. The table confirms what can also be inferred from 
the figure, namely that, for both tt and mesons /+ and nearly coincide. Comparison 
with the K form factors shows that the main SU(3)-breaking effect is in the normalization 
F(0), whereas the q^ dependence is only slightly modified. This can be understood from 
the fact that the formation of a vr or iiT meson is proportional to their respective decay 
constants fT^^Xi so that one would naively expect an enhancement ^ fx / fn = ^■'^ of the K 
form factors (at least if the three-parton states are not important), which is essentially what 
I find. 

Varying all input parameters within their respective allowed ranges, I obtain uncertain- 
ties between 5 and 10%. Combining this with the systematic uncertainty ~10% introduced 
by the need to separate the ground-state B contribution from that of higher states, the final 
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Figure 1: Form factors from light-cone sum rules in various approximations. 
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mm,) f^{m,) 


n PI /f 1 


13 


1 


m 

bp 


0.305 0.341 
1.29 1.41 
0.206 0.406 


^fm =/^(o) 

0.266 0.410 
-0.752 -0.361 


0.296 0.374 
1.28 1.42 
0.193 0.434 


0.27 ±0.05 0.33 ±0.05 
1.50 1.14 
0.52 0.05 



Table 1: Results for form factors with m;, = 4.8 GeV, sq = 33.5 GeV^ and = 6 GeV^ in 
the parametrization of Eq. (|12]). Renormalization scale for fp is n = rrib- The theoretical 
uncertainty is ~ 15 %. 





/;,latt(?')i i /T,LCSr(?') 


/o^latt(?^) 05 H /(^LCSr(?^) 


14.9 GeV^ 
17.2 GeV^ 
20.0 GeV^ 


0.85 ±0.20 0.85 ±0.15 
1.10 ±0.27 1.1 ±0.2 
1.72 ±0.50 1.6 


0.46 ±0.10 0.5 ±0.1 
0.49 ±0.10 0.55 ±0.15 
0.56 ±0.12 0.7 



Table 2: Comparison of lattice results for i? ^ tt form factors with results from light-cone 
sum rules. The errors for lattice results are those quoted in 0. 



uncertainty of the form factors is ~15%. A slight reduction of this uncertainty may be pos- 
sible if more accurate information on the twist 2 distribution amplitude becomes available, 
for instance from lattice simulations. 

A comparison with lattice results from the UKQCD collaboration is given in Table 2. 
The agreement with the lattice data is excellent, as it was also found for B ^ p, K* form 
factors in |jl6|. The LCSR point at = 20 GeV^ is just for illustration, because of which I 



also refrain from assigning it an error. 

5. In the limit m;, — * oo, Isgur and Wise have obtained a relation between the semileptonic 
and the penguin form factors [^: 



which is strictly valid only near zero recoil (i.e. near = ml). In Fig. 2 I plot the left- and 
right-hand sides of Eq. (|T3p over the full range of g^. The agreement between the curves 
for all is striking; they differ by 3% at = and by less than 1% at q^ = 17GeV^. A 
closer inspection of the underlying light-cone sum rules shows that to twist 2 accuracy 
is valid exactly and for arbitrary m, already at the level of correlation functions n± and 
Ft and thus — to that accuracy — is independent of the details of the extraction of the B 
meson contribution. For the twist 3 and 4 contributions to II-i- and Fp, (^) holds in the 
kinematical regime characteristic for the Feynman mechanism, i.e. near u ^ 1, and up to 
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Figure 2: Isgur-Wise relation (|T3D between penguin and semileptonic form factors. 



terms which are suppressed by one power of mb, which account for the 3% breaking of (|T3D 
by the hght-cone sum rule results. 

6. As the uncertainty associated with the twist 3 two-particle distribution amplitudes and 
the higher twist amplitudes is small, and the form factor thus depends essentially on the 
quality of information on the twist 2 distribution amplitude, it is worthwhile to investigate 
in more detail the impact of truncating the series @ after the first few terms. Let me first 
recall that in 03 (1 GeV) was determined from a QCD sum rule for Gegenbauer moments 
of in principle arbitrary degree, whereas aj was obtained from requiring 0,r(l/2, 1 GeV) = 
1.2 ± 0.3, which follows from a QCD sum rule for the vrA^A^-coupling. The reason for 
not considering QCD sum rules for higher moments is that they show a strong divergence 
with the degree n, rendering the series in (P) highly divergent. In order to simulate possible 
effects of higher moments without distorting the asymptotic u{l — u) behaviour near the end- 
points too much, I allow for a logarithmic divergence of the sum in Gegenbauer polynomials, 
yielding the following models: 

model I : 07^(1/2, /i = 1 GeV) = 1.2, model-dependent: 



^i{u, ^ = fib) = 6m(1 — u) — 0.95 • 6m(1 — u) - lnMln(l — u), 

5 



model II 



()^(l/2,/i = IGeV) = 1.2, 02(1 GeV) = 0.44 : 



(14) 



^l\u, /i = /ife) = 6m(1 -u){1 + 0.35 C2^\2u - 1)} 



■3/2/ 



-8.5-6u(l-n) |- lnnln(l-M) + l 

The corresponding Gegenbauer-spectra fall off as 

(I) : {a^} = {1,0.13,0.030,0.011,...} 



50 



Cf'(2u 



(15) 
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Figure 3: Dependence of /I on the twist 2 tt distribution amplitude. 



(II) : {a„} = {1, 0.35, 0.27, 0.10, 0.049, 0.027, . . .} 

In Fig. 3 I plot the above model distribution amplitudes as well as the one suggested by 
Braun and Filyanov p6|, i.e. 07r(4)- Although they look rather different, the resulting form 



factors also shown in Fig. 3, vary by at most 10%, i.e. are within the theoretical error. It 
is evident that the form factors do not depend on the details of the Gegenbauer-spectrum, 
but are sensitive only to a few gross characteristics like the value of the distribution am- 
plitude at one point (different from the end-points) and the first one or two moments. 
This is true as long as the distribution amplitudes are folded with smooth functions (as 
it is the case for form factors), so that higher order oscillatory Gegenbauer polynomials 
are effectively "washed out". A determination of the relevant few characteristics from an 
independent source, e.g. lattice simulations, would evidently help to further increase the 
accuracy of form factors calculated from light-cone sum rules. 

7. Summarizing, I have calculated the semileptonic and penguin form factors of B ^ it 
and B K transitions from light-cone sum rules. A new feature was the inclusion of one- 
loop radiative corrections to the leading twist contributions. The results are summarized 
in Fig. 1 and Table 1. The impact of radiative corrections and higher twist contributions is 
small, so that the achievable accuracy is limited by the inherent systematic uncertainty of 
light-cone sum rules, which is associated with the extraction of the B meson ground-state 
contribution out of the continuum of states coupling to the same current. This uncertainty 
is estimated to be ~ 10% and of the same size as the uncertainty induced by the input 
parameters in the sum rule. Hence, further refinement of the calculation by including 
higher twist contributions or two-loop radiative corrections is not expected to yield higher 
accuracy of the result. It would, however, be useful to have an independent determination of 
the few lowest moments of the twist 2 vr and K meson distribution amplitudes from lattice 
simulations. The existing results |2^ have large uncertainties, and in view of the recent 



improvements of the methods of lattice QCD and the availability of much more powerful 
computers, more accurate results seem within reach. Very recently |0 , a new method was 
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suggested to calculate the leading twist distribution amplitude on the lattice directly as a 
function of u. If feasible with statistical and systematic errors in the 20% range, this would 
help to reduce the total uncertainty of the B ^ n,K form factors to ~ 10%. The application 
of these lattice results would not be restricted to B meson decays, but also of direct relevance 
to the description of other hard exclusive processes, for instance single-meson production at 
HERA. 
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